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Radiative effects in the electroproduction of photons in polarized ep-scattering are calculated in 
the leading log approximation and analyzed numerically for kinematical conditions of current mea- 
surement at Jefferson Lab. Radiative corrections to the cross sections, their azimuthal distributions 
and Fourier coefficients are in particular focus. Kinematical regions where the radiative corrections 
are considerable are identified. 
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I. INTRODUCTION 

The processes of the photon electroproduction are in- 
tensively investigated both theoretically [l|, [2| and exper- 
imentally 3- 5] . The cross section of the process is sen- 
sitive to the deep virtual Compton scattering (DVCS) 
amplitude that is of great interest due to its connection 
to generalized parton distributions. The Bethe-Heitler 
(BH) process is not distinguishable from the DVCS mea- 
surements and therefore it is the basic background contri- 
bution to the observed cross section. One obstacle in the 
analysis of vast data on DVCS collected in Jlab experi- 
ments is the deficit of comprehensive theoretical calcula- 
tions of radiative corrections (RC) including the effects 
of hard photon emission with controlled accuracy. Avail- 
able calculations of QED radiative effects in [a-tlG] have 
certain limitations and cannot cover all modern require- 
ments of experimental data analysis on photon electro- 
production. In this paper we present the radiative cor- 
rection calculations to the cross section of BH in leading 
approximation. The main process contributing to the RC 
is two-photon emission, i.e., e+p — > e' +p' + 27. Another 
contribution is due to one-loop effects in e+p — > e'+p'+j. 
In the approximation the only leading term containing 
L = log(Q 2 /m 2 ) (m is the electron mass) is kept. For 
Jlab kinematics L ~ 15 and therefore, the used approxi- 
mation allows to keep the major part of the RC. 

The paper is organized as follows. The BH cross sec- 
tion is calculated in Section [TTJ Specific attention is paid 
on explicit representation of the BH cross section in- 
cluding polarization part of the cross section and mass 
corrections, as well as for angular structure of the BH 
cross section. RC calculation is performed in Section 
IIIII First we calculate the matrix element squared and 
trace all sources of occurrence of the electron mass de- 



pendence. Second, we represent the phase space of two 
final photons, introduce the so-called shifted kinemat- 
ics, and calculate integrals over additional photon phase 
space. Third, we add the contribution of loops and calcu- 
late the lowest order RC to the BH cross section. Fourth 
we generalize the result for the RC to the BH cross sec- 
tion to represent the higher order corrections. Section HVl 
presents the numeric estimates of the radiative effects in 
current experiments at JLab focusing on the RC to cross 
section in a wide kinematic region and angular structure 
(i.e., respective Fourier coefficients) with specific focus 
on the coefficient not appeared in the BH cross section 
but generated by RC. Finally in Section [V] we discuss the 
most interesting features of our findings in the theoret- 
ical calculation and numeric results (e.g., importance of 
mass corrections and kinematical regions with large ef- 
fects generated by RC), briefly describe the state-of-art 
in calculations of RC to exclusive photon electroproduc- 
tion processes and the place of our calculation among 
other calculations, and comment perspectives in further 
theoretical development in RC in exclusive photon elec- 
troproduction measurements. 



II. THE BH CROSS SECTIONS 



The BH process 

e(fci) +p(p) 



e'{k 2 ) +p'{p') +l(k), 



(1) 



is traditionally described using four kinematical vari- 



ables: Q 



-(fci - k 2 )\ x = Q 2 /(2p(k 1 - k 2 )), t 



(p — p') 2 , and <f), the angle between planes (ki,k 2 ) and 
(q,p') (<? = k x - k 2 ). 

The BH matrix element is Mbh = e 3 t" 1 J^ l J / f ff with 



' igor.akushevich@duke.edu 



J£ = u(p') ( 7^1 + Ur^-J^F* Up) (2) 



2M 



2 



and 



combinations of nucleon formfactors: 



u 2 



7m" 



fei — k ■ 



-U-2 



-2kki 

k\t k 2 e 
kk\ kk 2 



m k 2 + k + m 
e + e — 7 M 



7m 



2kk 2 

l^ke _ ekr/fj, 
2kk\ 2kk 2 



ui,(3) 



where u 2 = u(k 2 ), U\ = u(k±), and e is the photon polar- 
ization vector. The matrix element Mbh corresponds to 
the graphs in Figure [T^i and [I}). 

The cross section of the BH process is 

d^ = ^Ml H dr = 3 -^(j^ H ) 2 dr , (4) 

where S = 2ME\, E\ is the beam energy in the lab. 
system, and M is the proton mass. 

Phase space for the BH cross section is parametrized 

as 



1 d 3 k 2 d 3 p' d 3 k , 

dTo = Wf^W, 2^7 ( 1 

_ Q 2 dQ 2 dxdtd(j) 

f 2n 2 



p'-k) 



(5) 



with Ay = Si + AM 2 Q 2 kdAS x = S-X = Q 2 /x. Kine- 
matical limits on t are defined as 



-1 



((S X -Q 2 )(S X ± v / a7) + 2M 2 Q 2 ), (6) 



where W 2 = S x - Q 2 + M 2 . 

The 4-fold BH cross section (ctbh = d<r /dQ 2 dxdtd<f)) 
including both unpolarized and spin dependent parts is 



-o?Q 2 



AkS 2 xHJX^ 



^(T; + 2m 2 f i )7V 



(7) 



Those terms proportional to the lepton mass squared 
(m 2 ) are kept that give finite (i.e., non- vanishing for m — > 
0) contribution after integration over <f>. First two terms, 
{{T 1 + 2m 2 f 1 )F 1 and (T 2 + 2m 2 f 2 )F 2 ), describe unpolar- 
ized cross section and last two terms, ((T3 + 2m 2 T 3 ))F 3 
and (T4 + 2m 2 T4)J r 4), correspond to the spin dependent 
part of the cross section. The quantities T% are squared 





a) b) 
FIG. 1. Feynman graphs of BH process 



J r i=J r 3 = (F 1 (t)+F 2 (t)) 2 , 

^ = + 4^*?(«)). 

K = ~(Fi(t) + F 2 (t))(Fi(t) + ^F 2 (t))- 

The quantities are 
2 



(8) 



Ti = 
T 2 = 

T 3 = 



-{u 2 + w$-2Cyt), 



UqWo 

M 2 t 
-ttT 1 + 

2 u w 



(S 2 + X 2 



2 S X - Sw - Xu ), 
(2A(w - Q 2 ) - 2S(w + Q 2 ) 



AMt(rip') 
X t u w 

+ (UQ+W )(Q 2 - t)) 



-M 2 T 3 



+ —((Q 2 -uo)(t( V k 2 )+X(r ] p')) 

U W V 

(Q 2 +w )(t( v k 1 ) + S(r,p'))), (9) 



where X t = t(t - AM 2 ) and rj is the target polarization 
vector. The quantities represented the lepton mass cor- 
rections are 



f 2 
f 3 

T A 



2t 2t 



"0 

M 2 



s 2 + st x 2 + xt 



^MMl ( (2s + t)(— + — 

t-AM 2 \ {Zb+t) \u 2 + Sw 2 

+ 1 i(-Q 2 + } s(x 2 + {x-t) 2 )). 



M 2 f 3 + 2M(^(t( V k 2 )+X(r 1P >)) 
J_ + J_ )(t(^ 1 ) + 5(r ?P 0)). (10) 



All variables including the scalar products rjk\, r]k 2 , 
and r\p' are ultimately expressed in terms of 5 kinematical 
variables: S, t, Q 2 , x and <fi, e.g., 

w Q = 2kk x = -i(Q 2 + t) + ^(S X (Q 2 -t) + 2tQ 2 ) + 



A/; 



Ay 

u = 2kk 2 = w + Q 2 + t 



with S p = S + X and 



(11) 



Kw = AQ 2 W 2 (SX-M 2 Q 2 -rn 2 \ Y ){t-ti)(t 2 -t). (12) 

An additional azimuthal angle is required to describe the 
case of transversely polarized target (see eq. (fI3|l ). Note 
that in massless approximation (for m — > 0) the BH cross 



section exactly coincides with results of [ll| . The follow- 
ing equations relating our notation to the notation of ref. 
TjJ are valid: u = V 2 Q 2 , w = -ViQ 2 , and (for m -> 0) 
Kw = 4Q 4 S 2 S 2 K 2 . 

Explicit expressions for the scalar products of mo- 
menta with target polarization are 



r]ki = - 

V k 2 = - 

vp' = - 



SS X + 2M 2 Q 2 

2M-y/Ay ' 
XS X - 2M 2 Q 2 

-tS x + 2M 2 (Q 2 



t) 



2MJX 



(13) 



Y 



for longitudinal part of proton polarization vector (i.e., 
r] ||q) and 



rjki 
r\k 2 
VP' 



C0s(( / 9 + 0), 



As 



A' 



cos((p + ■ 



Q 2 SS X K 



: COS If 



(14) 



for its transverse part (i.e., r; 1 q). Here Xsx — 
SXQ 2 — M 2 Q A and <p is the angle between polarization 
and production planes, i.e., planes defined as (q, rj T ) and 
(q, p')). In the general case (77 = (0, r] x ,7]y, i] z )) the scalar 
products are 



i-jki = 
r]k 2 = 
VP = 



sx 



SS X + 2M 2 Q 2 

XS X - 2M 2 Q 2 
2MV\T 



\5X 



Vx 



Q 2 SS X K 

V^SX^Y 

-tS x + 2M 2 (Q 2 - t) 



(r] x cos (fi + rjy sin t 

Vz- 



2MV\^ 



(15) 



The cross section <jbh is defined as the 4-dimensional 
cross section in eq. Q • It means that integration over tp 
is assumed to be performed resulting in additional factor 
27r. This is because of respective symmetry of the process 
with unpolarized or longitudinally polarized target. The 
cross section of transversely polarized target explicitly 
depends on tp because of (fl4| . Thus u in the case of 
transversely polarized target means 5-dimensional cross 
section with addition of dtp/2ir. 



A. Angular structure of the BH cross section 

Azimuthal structure of the BH cross sections are often 
of interest both theoretically and experimentally. Depen- 
dence of the BH cross section (|U) on the angle </> appears 
in uq and wq in numerator and denominator of the BH 






a) b) c) 

FIG. 2. Feynman graphs of two real photon emission 

cross section and in scalar products r\k\ and r]ki for trans- 
versely polarized proton (angle ip is assumed to be fixed, 
i.e., ^-independent). The upolarized and spin-dependent 
parts of BH cross section can be presented in the form: 

/ 



unp 

T BH 



r LP 
'BH 



TP 
'BH 



f 

V1V2 
f 

V1V2 



cos 20), 

{COXP + Ci,lpCOS0), 



cq,tp + c\.tp coscj) + si.tp sin0), (16) 



where / = a 3 S , 3/(87ra; 3 a^ / ' 
are expressed as 

1 



The Fourier coefficients 



c " = 7T~F / d(j)ViV 2 cr B H, 
2nj Jo 



Cl 



1 

t/ Jo 
1 



)COS0 V1V2 CT B H, 



2tt 



c 2 = — / dcj) cos 2(f) V1V2 <Jbh, 
Jo 



2tt 



si = — / d</)sin0 V\V 2 <Jbh- 
Ti"/ Jo 



(17) 



Only terms T\ — T4 contribute to eq. (fT7| while mass cor- 
rections represented by T± — T4 can be neglected. This is 
because the integration over <fi in pT|) is performed with 
weights V\V 2 reducing singularity level in this terms re- 
sulting in their zeroth contribution to (I17[) in massless 
approximation. The Fourier coefficients (|17[) are defined 
in exactly the same way as those given in eqs. (35-42) of 
11]. Our analytic calculation of the Fourier coefficients 
for unpolarizad, longitudinally and transversely polarized 
cross sections and their subsequent analytical compari- 
son with the expressions of |11J show that both sets of 
formulae are identical. Therefore, we do not show the ex- 
plicit expressions for the Fourier coefficients here. Note, 
however, that for transversely polarized case the Fourier 
coefficients still depend on ip (defined after eq. [T4"|) that 
must be assumed to be fixed to have eqs. (fTB"]) and (TP7| 
valid. Alternatively, one can assume that the angle be- 
tween scattering and polarization planes {(p. = tf> + ip). 
In this case ^-dependence for transversely polarized tar- 
get needs to be moved out from the expressions for the 
Fourier coefficients in Eq. (fT6")l and that equation needs 
to be rewritten as 



TP 

a BH 



f 



■ (Cg tp COS (fi + c[ TP COS ip COS ( 



(18) 



4 



/ 

+ s[ T P sin ip sin rp) = — ((c 1TP - s[ Tp ) cos tp 
l F\ F 2 

+2c' 0TP cos(</> — <p) + {c' 1TP + s'i.tp) cos(20 — p j) 
f 



V X V, 



■ (c tp COS tp + c 1TP cos(0 — p) 



-c 2 tp cos(20 — tp)). 



Important research questions are what magnitude of 
RC for Fourier coefficients is and whether RC can gen- 
erate new functions (e.g., involving cos 30 or sin 20) van- 
ishing at the level of the BH cross section. 



III. RC CROSS SECTION 

The cross section of two photon emission, i.e., the pro- 
cess 

e(fci) + pip) — > e'(*a) + p'{p') + + 7(^2), (19) 



is 



da 



4S* 



!=1 



(20) 



dT = 



where additional factor 2 in the denominator is because 
there are two identical particles (photons) in the final 
state. Phase space is parametrized as: 

1 d 3 k 2 d 3 p' d 3 Ki d 3 K 2 
(2?r) 8 2E 2 ~2p~^ 2wi 2uj 2 
x8 i {k x +p-k 2 -p l -K1-K2). (21) 

Six matrix elements of the process with emission of ad- 
ditional photon correspondent to graphs in Figure [5] are 
denoted Mi-e = e 4 i _1 Ji_6,ju- The quntities </i-6,^, 
proportional to the leptonic currents, are: 



J 



fci — k + m „ k\ — k\ 



■2nki + V 2 



C-2- 



-2k\K\ 



-eiui, 



2 M = "27m 



fci — k + m „ fci — k 2 + m , 



ei- 



= u 2£2 



-2nk x + V 2 -2k lK2 
k 2 + k 2 + m „ k 2 + k + m 



-£2"i, 



2k 2 K 2 



-<"i 



2 K fc 2 + V : 



-7m u i> 



J. 



4 M = "2£l- 



&2 + ^1 + m . fc 2 + k + to 



2fc 2 Kl 



-£2 



2 K fc 2 + y ; 



-7m U 1' 



5^ = U2£l- 



ki + m k\ — k 2 + to „ 



2fc 2 Kl 



— 2k\K 2 



-e 2 u l7 



_ „ k 2 + k 2 + m k\ - ki + to , 
J6p = "2^2 ^777^ In —, eiwi, 



2fc 2 K 2 

where V 2 = n 2 = (k\ + n 2 ) 2 



-2k\K\ 



(22) 



A. Matrix elements in leading approximation 

There are four kinematical regions contributed to the 
cross section in leading approximation: when one of the 



photon is observed and another in so-called s- and p- 
peaks. For s-peak (p-peak) the additional unobserved 
photon is emitted in the direction of the initial (final) 
lepton. Therefore, 



M 2 



2 s 



(23) 



where indices correspond to the unobserved photon, e.g., 
Is means that the photon with momentum k\ is unob- 
served and in the s-peak. 

The matrix element A4 2 s squared in the leading ap- 
proximation is calculated assuming that the momentum 
Ki of unobserved photon is approximated as 



Hi = (1 - Zi)fci. 



(24) 



However this approximation have to be carefully applied 
after analyzing the the structure of poles, i.e., powers 
of fci«i in denominators. Only terms with the first-order 
pole (1/fciKi) contribute the the cross section in the lead- 
ing approximation. The second-order poles appear in the 
form of m 2 / (kiKi) 2 and does not contain the leading log 
after integration and taking the limit to — > 0. Only J Xfl 
and Jq^ have the pole, 



J 



Jt 



6/i 



2 k\K 2 k\K\ 

Z\ k\t\ 
2 k\K 2 k\K\ 



U2l^{z\ki - k 2 )e 2 ui 
u 2 e 2 (k 2 + k 2 )j^ui, 



(25) 



while J 2m , J3 M , J4 M , and do not. It means that the 
interference term (J lM + Je^){J 2u + <hv + J<kv + Jn^Y has 
the pole and therefore contributes to leading log approx- 
imation and that + Jg M squared can have the pole of 
the second order. 

Calculating the interference, the equations (|2"5]) can be 
applied, and the substitutions (IM|) and m — can be 
used everywhere except in k\K,\ in the denominator. Both 
Ji M and Jq^ are proportional to kiei. Therefore, to cal- 
culate the interference one needs to calculate Ji^\^\ 
(i = 2, . . . , 5) by averaging over unobserved photon po- 
larization states. This results in 



u 2 e 2 {k 2 + k 2 )-i^ui 



2(1 - z\) k 2 K 2 



E, / j j \ w 2 7 M (ziA;i - k 2 )e 2 ui 
kiei(J 2fl + J 5m - — — — . 
2 k\K 2 zi(l - zi) 



(26) 



and therefore 



(Jifj, + •h i i){J2v + -hu + Jiu + -hv)^ + h.c. = 

= J^ H {z 1 k u k 2 ){J? H {z 1 k 1 ,k 2 ))^- h ^ V ( 2? ) 

felKl (1 - Zi) 



Calculating the term with Jxu + J&n squared, the 
poles have to be extracted in the form of \fk\K,\ and 
TO 2 /(fci«i) 2 , and only then the substitutions (|2~4"|) and 



5 



<3 Z 

e z 4 


p' 










e'/ / 

fly , 








y 



FIG. 3. Definitions of vectors and angles in the Lab. frame 



m = can be used everywhere except in k\K\ in the 
denominator. This results in 

(Jlfi + J&^iJlu + ^6f) t Ui^(l-zi)fei 

= J* H ( Zl h, k 2 )(J^ H ( Zl k U fc 2 ))t- 1 ' 1 



Z\ K\K\ 

We finally have in leading approximation: 



M{ 



1 AmOL ka2 



km_ 



M BH (ziki,k 2 )- 



1 + zf 



(28) 



(29) 



zi(l - Zl) 

In the case when the unobserved photon is emitted 
parallel to the final electron the scalar product fc 2 Ki is 
small. For this case it is assumed that k\ = (z^ 1 — l)k 2 
resulting in 



V Z 2 J 1 - Z 2 



(30) 



B. Phase space and shifted kinematics 

The photon four-vectors appear in denominators of 
(|29|) and (|30[) in the form of scalar products kiKi and 
k 2 ni. Two integrals over phase space of two photons are: 

d 3 Ki d 3 K2 S(A — K\ — K 2 ) 

2uj\ 2oj 2 k\K\ 

d 3 Ki S(A 2 - 2A/ti) _ ttL 
2loi k\K\ w ' 

d 3 Ki d 3 K 2 5(A — Ki — K 2 ) 

2cji 2w 2 fc 2 Ki 

d 3 Ki 6(A 2 - 2Aki) _ ttL 

2l0\ k2Kl u ' 

where A = k\ + p — fc 2 — p' , w = 2k\ A, u ~ 2fc 2 A. Only 
terms containing the large (or leading) logarithm L are 
kept. The results (|3"Tj) are immediately obtained if to con- 
sider the system of center-of-mass of two photons (A = 0) 
with z-axis directed along ki and k 2 , respectively. 



The phase space of the final proton is parametrized as: 
d 3 p' 



8M 



A* , , , a , dtd^dV 2 
— dtdcpd cos 9 = 



4aA 



Y 



(32) 



where V 2 — A 2 and 9' is the angle between q and p' 
The relation 



V 



tS x 



+ t- 



■ COS C 



(33) 



2M 2 ^ 2M 2 

was used to obtained the parametrization in terms of V 2 . 
Finally the integration over dT is 

dT „ L 



dT 

k 2 Ki 



dT 



dT n 



8n 2 w 
ttL 



dV 1 



-dV z 



(34) 



The matrix elements squared for s- (and p-) peak con- 
tributions in eqs. (|29|) and (|30]) are expressed in terms of 
zi and z 2 , therefore the variable V 2 (and cos 9') has to 
be related to these variables. The equation for establish- 
ing this relation is obtained from condition in ^-function 
argument of intermediate expressions in (|31[) if to use 
the representation for n\ used in subsection IIII Al i.e., 
Ki = (1 — zi)ki for s-peak and K\ = (z^ 1 — l)k 2 for p- 
peak. Below for representation of this equation and its 
solution we use the generalized notation included both 
z\ and z 2 . Substitution z 2 = 1 [z\ = 1) has to be used 
to formally extract s-peak (p-peak) contribution. Also 
we define the 4- vector q z : q z = z\k\ — fc 2 for s-peak, 



ki 



_1 fc 2 for p-peak, or q z 



Zlki 



z 2 ~ l k2 m 



the generalized notation. Meaning of the used vectors is 
clarified in Figure Vector q z has meaning of "true" 
transferring momentum in the case of additional photon 
emitted. The vector is in the plane OXZ, its projection 
into OX and OZ axes are always negative and positive 
respectively. Its magnitude is always less than that of q. 
The equation for establishing the relation between z\, z 2 , 
and V 2 in terms of introduced notation reads: 



A 2 - 2A{q - q z ) = 



Yz 



2M 2 



(cos0-A) = O, (35) 



where cos 9 is the angle between q 2 and p'. It is expressed 
in terms of the angle between q and q z (denoted by 9 Z ) 
as 



cos tt = cos f cos f 



smf sintf* cos< 



(36) 



where sinus and cosine of 9 Z and the quantity A are de- 
fined by kinematics in terms of z\. 2 and measured quan- 
tities: 



cosf 2 = 



sin u z = 



S x { Zl S - z^X) + 2(z~ 1 + Zl )M 2 Q 2 
2(z- 1 - Zl )MQ(SX - APQ 2 ) 1 ' 2 

\/Ay\/ Xyz 

(ziS - z~ x X)t + 2M 2 (t - z^Q 2 ) 



A = - 



\ Yz = { Zl S - Z2 Xy + AM 2 Zl z 2 




(37) 



Eq. (|35[) has unique solution in the kinematically allowed 
region: 

Acos8 z + \/^osin( 



cos 9' = 



COS( 



sin 9' 



T> 



cos^ 
cos 9 Z 



9 Z + sin 2 Z cos 2 
FDq — A sin 9 Z cos ( 



2 z + sin 2 6> z cos 2 . 



cos 



cos 

2 



sm 



: cos 



A 



(38) 



The direction of q z defines new polar (9) and azimuthal 
((/>) angles of the final proton and thus generates so-called 
shifted kinematics. The angle (9) was defined in (|36|) and 
the angle </> is related to measured <f> as 



cos V, sm V cos ( 



sm o, cos ( 



and 



sin <b sin 9 = sin 0' sin < 



(39) 



(40) 



The origin of the equation (|4T)|) is clear because the pro- 
jection of p' on axis OY is the same for original and 
shifted kinematics. Recall, that the equations (|35M0j) 
are used for both s- and p-peaks. In the first case one 
sets Z2 = 1 and Z\ = 1 is set for the second case. 

The target polarization in shifted kinematics is calcu- 
lated using the orthogonal transformation: 



Vv 

Vz 



cos9 z i] x + sm9 z r/ z , 



Vy, 

— smt 



z rj x + cos 9 ' z rj z 



(41) 



Thus, in the shifted kinematics the target polarization 
are not longer pure longitudinal or transversely polarized, 
therefore, the scalar products in shifted kinematics are 
then calculated using eqs. (p~5]) . 

The variable V 2 is related to Z\p, through eq. (|33|) 
where COS0' is given by (|3"5)l and quantities in the R.H.S. 
of (138)) depend on zx,2 in (|37)) . Tedious, but straightfor- 
ward calculation gives 

1 dV 2 ~ 



fXy sin 9' 



w dzi 
IdV 2 




dz 2 



(42) 



Z 2V A Yz 

The R.H.S. of the equations (I40I42|) are taken using re- 
spective peak kinematics. 

The equation for minimal value for Z\ (denoted by z™) 
allowed by kinematics is cos9 z = A. It follows from 
(33): V 2 ax = tS x /2M 2 + t-Q 2 + y/)^y/X~ t /2M 2 (that 
corresponds to cos (9' = 1). The solution is 



Xt - 2M 2 t + £(XS X - 2M 2 Q 2 ) 



where £ 2 = 
we obtain 



St- 

At/Ay 



(43) 



2M 2 Q 2 + £[SS X + 2M 2 Q 2 ) ' 
Similarly, for kinematics of p-peak 



m = Xt + 2M 2 Q 2 + £(XS X - 2M 2 Q 2 ) 
7 " 2 ~ St + 2M 2 t + i(SS x + 2M 2 Q 2 ) ' 

Both zY 1 and zl 71 do not depend on <fi. 



(44) 




FIG. 4. Feynman graphs of one-loop effects for the BH cross 
section 



The lowest order RC to BH cross section 



Combining results obtained in Sections IIII Al and IIII Bl 
we find the cross section of two photon emission as: 

a s (S,x,Q 2 ,t,(t>) = ^-L x 



2tt 



dzi 



1 + z\ sin 



1 - z i V, 



\jj (— ) ctbh(ziS, x s ,ziQ 2 ,t, (j) s ), 



O.s 



<Tp(s, x, Q 2 , t, 



2tt 



-L x 



(45) 



dz 2 



l + z? sin 



2fe(l - z 2 ) V, 



T/jy-^r) <Jbh(S,x P7 z^ Q 2 ,t,(j) p ), 



Op 



where x s = z\Q 2 j(z\S — X) and x p = Q 2 j{z-iS — X) are 
Bjorken x in shifted kinematics; sin#' and <f) are given by 
(|3"5|) and (|4"0")) — subscript explicitly indicates the type of 
kinematics for that these quantities have to be calculated. 

The integrals in (|45[) are divergent at upper integra- 
tion limit, therefore it is regularized using a parameter 
w m in separating the integration region on the part corre- 
sponding to emission of soft and hard photons. For zx,2 
the regulating parameter A is A — Ai = 2Mui m i n / S for 
s-peak and A = A2 = 2Mw m i n j X for p-peak. 

The contributions of loops (Fig. HJa-Zi) and soft photon 
emission are known Q . Their sum is proportional to BH 
cross section 



aV =n{ l ° g ^X— ~ 
and can be presented as 

l-Ai 

aL ( f , l + z 2 
cry = - — <tbh\ / dZx- 



Lubh 



(46) 



i-a 2 



2tt 



1-zi 



1 + 4 

1-Z2 



(47) 
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Sum of (145)) and (|47|) is infrared free and regularization 



can be removed: Ai 2 
cross section is: 



0. The result for the observed 



^obs°° P (S, x, Q\ t, <f>) = (1 + 2Il(t))a BH (S, x, Q 2 , t, 0) + ^-L 



1 + zj' 

1-Z2, 



sin 0L 



1/2 

Op 



2tt 




?r) (^r) obh{ziS, x Sl ziQ 2 , t, 4> s ) - ct B h{S, x, Q 2 ,t, <p) 

1 / X \ ^ 1 — 

J 2 1 )— (— ) <Jbh(S,x p ,z2 Q 2 ,t,(f)p) - <7 B H(S,X,Q 2 ,t,< 



(48) 



Here II(i) = a/(2ir)S vac and <5„ ac is the contribution of 
vacuum polarization by leptons and hadrons (Fig. 
calculated as in [l2j (see eq. (21) and discussion before 
cq. (20)). 

Behavior of the cross section for t close to kinematical 
bounds (i.e., in the region where t ~ t\ and t ~ tz) 
deserves special attention. The integrals in (|4"5|) become 
infinite when t — > t\ or t — > t 2 . In this limit z™ = 1 and 
z™ = 1 . To extract the divergence the part of integrals in 
(US]) from to z™ or z™ need to be calculated analytically 
resulting in: 



(Sva 



Jinf 



Ifin 



VBH + &F- (49) 



where op is the non-divergent contributions of remaining 
integrals (i.e., as in (|4"5]l . but with low limits z™ and z™). 
The correction terms 



5 fm = j (z?{2 + zH + z™(2 + z™)) , 

Sinf = i(l0g(l - Zl) +l0g(l - Z 2 )). 



(50) 



represent the the finite and infinite parts of the results 
of the analytical integration. The source of occurrence of 
the divergence is known 13]. The divergence is canceled 
by taking into account multiple soft photon emission. We 
follow the so-called exponentiation procedure suggested 
in 0: 

(l + -(&uac + <W + £/in)) -> (51) 

exp(^<5 in/ ) (l + -(<W + <5/m))- 

After this procedure the observed cross section vanishes 
at the kinematical bounds on t. 



D. Higher order corrections 

In previous section we found RC to BH cross section in 
leading approximation induced by lepton leg in the lowest 



I 

order over a. The generalization of eq. (|48[) on highest 
order over a using electron structure function method 
suggested in [l5| (see also [H, GjJ) has a form: 



°"ofcs {S, x,Q , i, 



1 1 

= J dz x J ^lD( Zl ,Q 2 )D(z 2 ,Q 2 ) 



sin 



— U2°bh(ziS, x sp ,z 2 ziQ ,t,(f>), (52) 

x l X>q' 

where a BH = vbh[o? -> a 3 /(l - n(i)) 21 
ziQ 2 1 \z\Z2S — X) and 



( ■ Sp 



6 2,1 



Xt + 2 Zl M 2 Q 2 + £(XS X - 2M 2 Q 2 ) 
Zl St + 2MH + z^(SS x + 2M 2 Q 2 ) 



(53) 



The electron structure function D{z, L) includes con- 
tributions due to photon emission and pair production 



D = 



n e+e 

N 



D% +e ~ 



(54) 



where is responsible for the photons radiation and 

D N andD s describe pair production in non-singlet 
(by single photon mechanism) and singlet (by double 
photon mechanism) channels, respectively. The explicit 
expression for D(z, L) are given by eqs. (5-7) of ref. [l6| . 

Notice, that equation (|48p can be reproduced by ex- 
pansion of (1521) over a and keeping only zero and first 
order. 



IV. NUMERICAL ESTIMATES 

Numerical analysis is designed to evaluate the RC for 
the cross section and the Fourier coefficients in the kine- 
matics of modern measurements at Jlab (D-d. Specific fo- 
cus in this analysis will be on i) the t- and ^-dependencies 
of the magnitude of RC factor and ii) investigation of RC 
for the Fourier coefficients both non-vanishing and van- 
ishing at the level of the BH cross section. 
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RC factor 



3.5 4 , 
-t, GeV 



FIG. 5. The observed cross sections of the BH process (upper 
plot) and respective RC factors (lower plot) for beam energy 
5.77 GeV, x=0.4, and Q 2 = 1.8GeV 2 . The line with (without) 
dots shows the results of calculation without (with) the cut 
on missing energy (E~, <0.3 GeV). 



A. Cross section 

The i-distribution of the BH cross section has two 
sharp peaks that correspond to collinear radiation. The 
typical shapes of the ^-dependence of the BH cross sec- 
tion with RC are represented in Figure [5^,. Figure [SJd 
gives ^-dependence of the RC factor for the given kine- 
matical points. The plots for spin dependent parts looks 
similar for both longitudinal and transverse polarizations 
(not shown) . In this analyses the cross section integrated 
over <f> is considered. 

Analyses of i-dependence presented in Figure [5] re- 
vealed three specific regions in which the shapes of RC 
deserve attention and further clarification: i) the region 
close to bound over t where RC factor goes rapidly down, 
ii) the region close to collinear peaks, and iii) the region 
between the peaks where RC factor can reach large val- 
ues, however, capable of being suppressed by a cut on 
missing energy (also shown in Figure [5] by the line with- 
out dots). 

Decrease of the RC factor in the region close to the 
bounds (i.e., t ~ t\ or t <~ t 2 ) is simply the reflection of 
the fact that observed cross section after the exponentia- 
tion procedure (j5"Tj) . as well as the observed cross section 
(15^1) included higher order corrections, goes to zero at 
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t = -0.5 GeV 2 


t = -1.2 GeV 2 
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t= -1.5 GeV 2 


, t = -1.8 GeV 2 


Values of RC factor 

for <|> h = tt: are 
30 for solid line and 
36 for dashed lines. 





0.2 0.4 0.6 0.8 
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<s>J2% 



FIG. 6. The RC factor vs. r/> for beam energy 5.77 GeV, 
£=0.4, and Q 2 =1.8GeV 2 . The solid (dash) lines correspond 
the result of calculation without (with) the cut on missing 
energy (£ 7 <0.3 GeV). 



these kinematical bounds. 

In the region close to s- and p-peaks, i.e., when 

Q 2 X x ± Q 2 S 



t = U = - 



S-Q 2 



t / i. 



X + Q 2 



(55) 



the RC factor slowly decrease (when —t approaches to the 
s-peak from the left or p-peak from the right, see Figure 
[S]), reach its minimum at t — t s or t = t p , and then 
rapidly increase, reach its maximum at t = — Q 2 . Anal- 
ysis of the integrand shows that the region around the 
point <f) = 7r is responsible for this difference. Therefore 
the (^-dependence of the RC factor (Figure O was ana- 
lyzed. The RC factor typically has flat behavior except 
the point <p = n corresponding to the situation when the 
scattering and production planes coincide. In this case 
the RC factor can rapidly increase. Further analysis of 
the integrand showed that this increase of the RC factor 
is due to contribution of the second integral in (l4"5l) when 
Wo is very small. The the second integral in (|4"8"j) describes 
the p-peak contribution of the one photon, and the re- 
gion of small wo corresponds to the s-peak of the second 
photon. Therefore, the large contribution comes from the 
two photon emission process when two irradiated photons 
are collinear to initial and final electrons. Correspond- 
ing BH process (i.e., one photon emission process) is the 
process with the emitted photon with 4-momentum cor- 
responding to the sum of momenta of the two collinear 
photons. This photon is not collinear and therefore re- 
spective cross section of BH process is not large. The RC 
factor defined as the ratio of observed cross section (with 
large contribution of the two collinear photons) to the 
BH cross section (with not large BH cross section) can 
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become larger than 2, i.e., RC to BH cross section can 
be larger than the BH cross section. Roughly the effect 
for RC factor can be estimated as 1 + aL 2 (one collinear 
photon produce one leading log L). If L ~ 15 then the 
RC factor equals 2.64. 



B. RC and azimuthal structure of the cross section 

Azimuthal structure for the unpolarized BH cross sec- 
tion and for the longitudinally and transversely polarized 
cross sections are represented by Fourier coefficients de- 
fined in Section III Al There are eight non-zero Fourier 
coefficients: three for unpolarized cross section, two for 
longitudinally polarized, and three for transversely polar- 
ized. Radiatively corrected Fourier coefficients are calcu- 
lated using eqs. (|17I18I) with the observed cross section 
(1521 substituted instead of (Tbh- Figure [7] presents the 
results for these coefficients calculated using the BH and 
observed cross sections. The observed cross section was 
calculated with and without kinematical cut on maximal 
photon energy E 1 = 0.3GeV". One can see from this plot 
that the eight coefficients are quite stable in respect to 
RC. Similarly to the case of the cross section the regions 
with noticeable effect from RC are the region of small 
—t and the region of t close to (and between of) the s- 
and p-peaks. Also the results show that using the cut 
on missing energy suppresses the correction in the latter 
region. 

In contrast to the BH cross section, the azimuthal 
structure of the observed cross cannot be represented 
neither in terms of this eight coefficients nor in terms 
of any finite number of such coefficients. This is because 
of complicated and nonlinear dependence of the observed 
cross section on <f>. Several coefficients representing next 
terms in the Fourier series are presented in Figure [5] All 
of them are defined through cos(n(f>). The Fourier coef- 
ficients with sin(n</>) were also investigated (all of them 
vanish at the level of the BH process), however no signif- 
icant contributions at the level of observed cross section 
were found. 



V. DISCUSSION AND CONCLUSION 

In this paper we calculated RC to the BH cross section 
in leading approximation. Both unpolarized and polar- 
ized parts of the cross sections were considered. All final 
formulae are presented in analytical form. Details of cal- 
culation of matrix element squared are given with specific 
attention to occurrence of mass terms non-vanishing in 
the approximation of small lepton mass. Phase space 
was parametrized using the notion of shifted kinematics 
resulting in compact and convenient parametrization for 
the two-photon phase space and opportunities for ana- 
lytic integration over angles. Numerical analysis of the 
effects of RC was focused on the RC to cross section and 
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FIG. 7. The Fourier coefficient with nonzero contribution 
to the BH cross section (|16|) . Points represent the BH cross 
section, and solid (dashed) line represents the observed cross 
section (|52|l calculated without (with) the cut E 1 = 0.3GeV . 
The coefficients were rescaled to use the same scale on all plots 
and to compare them among each other: the original values 
of coefficients are calculated by dividing the values obtained 
from plot by the rescaled factor. 



the Fourier coefficients representing the angular depen- 
dence of the BH and observed cross sections. 

Analysis of the RC to the BH cross section revealed the 
kinematical regions where RC can exceed the BH cross 
section in several times. This is the region with scattering 
azimuthal angle equaling tt. The situation when both 
photons are collinear (one is collinear to initial lepton 
and another is collinear to final lepton) are kinematically 
allowable. Since the photon in respective BH process is 
not collinear (its momentum is the sum of two collinear 
photons), the BH cross section is not so large. As a result, 
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FIG. 8. The Fourier coefficients not contributing to the BH 
cross section (|16|) . Solid (dashed) line represents the observed 
cross section (|52[) calculated without (with) the cut E-, = 
0.3GeV. The definition of the rescaled factors are the same 
as in Figure [7] 



the RC factor can be at the level of several dozens. 

There are eight Fourier coefficients contributing to the 
BH cross section with arbitrary polarized target, while 
new coefficients appear in the Fourier expansion at the 
level of the observed cross sections. The calculation of 
the contributions of the additional terms of the Fourier 
series can be significant. For example, as follows from 
the comparison of results in Figures [7] and |8] the effect 
of "new" coefficient C3 yUnp can reach 10% from the ef- 
fect of main contribution represented by Co jUn p. This 
effect, however, can be suppressed by using the kinemat- 
ical cut on missing energy. Note that experimental pro- 
cedure of extraction of the Fourier coefficients from data 
is based on the fitting of observed cross section by the 
functions representing the angular structure of the BH 
cross section. Occurence of large terms of the next or- 
ders in Fuorier expansion can results in systematical un- 
certainties in kinematic regions where the effects of these 
additional terms is noticable. 

One feature of the calculation is that the lepton mass 
cannot be completely eliminated in the expressions for 
the BH cross section First, the lepton mass has to 
be kept in the lepton propagators wq and uq. Since the 
propagators are proportional to El — plcos/3 (El and 
Pl are energy and momentum of the lepton, and /3 is 
the angle between momenta of the lepton and photon), 
there are kinematical points where wq or uq vanish in 
massless approximation making the BH cross section in- 
finite. These points can be excluded when the BH process 
is investigated experimentally. However RC calculation 



requires integration of the BH cross section over broad 
kinematical area and the singular point occur in the in- 
tegration region. Therefore the lepton mass has to be 
kept in the expressions for wq and uo. This is the rea- 
son of occurrence of the m 2 in (fl2l) . Second, terms in 
the BH cross section containing m 2 in numerator and wq 
or uq squared in denominator are also infinite in massless 
approximation for certain <fi and result in finite (indepen- 
dent of the lepton mass) terms after integration over <fi. 
Such terms were kept in the expression for the BH cross 
section Note that our experience of dealing with RC 
tells us that such terms can give important contribution 
to the observed cross section (e.g., in DIS cross section 
measurements). 

The motivation for our calculation was the lack of com- 
plete calculations of the RC performed with accuracy to 
be controlled. One-loop correction and soft photon emis- 
sion was calculated by Vanderhaeghen et al. 0- Detailed 
consideration of one-loop correction was done. Box-type 
diagrams were evaluated in the style of ref. [HI . How- 
ever, the radiative tail corresponding to photon emission 
processes was calculated in the approximation where the 
photon energy is very small compared to the lepton mo- 
menta. 

Bytev, Kuraev, and Tomasi-Gustafsson @ applied the 
method of the electron structure functions to calculate 
RC due to two photon emission in the process e _ /x + that 
was chosen as a model process of DVCS. Main focus in 
this calculation was on the correction to the helicity-odd 
part of cross section, i.e., the interference between BH 
and DVCS amplitudes. Authors were focused on an- 
other experimental design: they integrated over the en- 
ergy fraction of scattered electron. 

The task of the calculation of RC to BH is closely re- 
lated to the task of RC to radiative tail from the elastic 
peak that is the important (and often dominant) contri- 
bution to RC in DIS measurements. The radiative tail 
is simply the BH cross section integrated over photonic 
variables, i.e., over <f> and t. Integration over <j) is per- 
formed analytically, and because of dependence of the 
cross section on formfactors the integration over t is left 
for numerical analysis. Programs for RC calculation of 
the radiative tail such as POLRAD 2.0 [l9| include both 
the contribution of the radiative tail (correspondent to 
the BH cross section) and approximate calculation of the 
RC to the radiative tail (correspondent to two-photon 
emission and loop effects) 0, [2(| ■ The approach to exact 
calculation of the RC to the radiative tail was developed 
by Akhundov, Bardin, and Shumeiko @. They used the 
formalism of covariant extraction and cancellation of in- 
frared divergence and calculated the QED corrections to 
the elastic radiative tail for unpolarized case. No an- 
alytical expressions represented the result of the exact 
calculation were published. 

The formulae in this paper are presented in analytical 
form providing good starting point for more precise cal- 
culations. One further generalization can be done using 
the approach of @ to exactly calculate the lowest order 
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correction to polarized BH cross section. Another direc- 
tion for generalization is to apply the developed formal- 
ism for RC to DVCS, i.e., interference of BH and DVCS 
amplitudes. Hadronic part of DVCS is known from refs. 

Urn 

Note also that the developed formulae are obtained 
for the specific way of reconstruction of kinematic vari- 
ables. Specifically, leptonic and hadronic momenta are 
used to reconstruct the kinematics of the BH process. 
Kinematical variables of the photon were assumed to be 
unmeasured. If information about photonic variables are 
involved into reconstruction of the kinematics of the BH 



process the calculation presented in this paper requires 
modification. Universal way to avoid multiple calculation 
to cover all possibilities for data analysis designs is the 
development of the Monte Carlo generator of the BH pro- 
cess with the additional process with two photons. Any 
specific choice of base set of kinematical variables can be 
used for this construction including those considered in 
this paper. 
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